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We theoretically investigate the dynamics of a trapped ion immersed in a spatially localized buffer 
gas. For a homogeneous buffer gas, the ion’s energy distribution reaches a stable equilibrium only if 
the mass of the buffer gas atoms is below a critical value. This limitation can be overcome by using 
multipole traps in combination with a spatially confined buffer gas. Using a generalized model for 
elastic collisions of the ion with the buffer gas atoms, the ion’s energy distribution is numerically 
determined for arbitrary buffer gas distributions and trap parameters. Three regimes characterized 
by the respective analytic form of the ion’s equilibrium energy distribution are found. Final ion 
temperatures down to the millikelvin regime can be achieved by adiabatically decreasing the spatial 
extension of the buffer gas and the effective ion trap depth (forced sympathetic cooling). 

PACS numbers: 


The ion motion inside a radio frequency (RF) trap 
is characterized by the interplay between a fast oscilla¬ 
tion driven by the RF-held (micromotion) and a much 
slower oscillation in the conhning ponderomotive poten¬ 
tial (macromotion) [T], thus representing a prototypi¬ 
cal example of a dynamically driven nonlinear system 
[5]. Through elastic collisions with a cold buffer gas, 
either consisting of a cryogenic noble gas [SHi] or laser 
cooled atoms ISHU], the ion’s motion can be efficiently 
reduced, thus opening a wide range of applications rang¬ 
ing from precision spectroscopy m and spectrometry 
[Mini over quantum computation [M to cold chemi¬ 
cal reactions dynamics [M and astro-chemistry [Mill]- 
However, elastic collisions influence the permanent ex¬ 
change of energy between micromotion and macromotion 
resulting in a net energy transfer from the micromotion 
to the macromotion [MIM- ^ consequence of this 
coupling, the ion’s hnal mean energy generally exceeds 
the buffer gas temperature [3 IHHM and the ion’s en¬ 
ergy distribution is predicted to deviate from a thermal 
distribution [M [Ml ES] • Due to the stability constraints 
of a RF ion trap, efficient cooling through collisions with 
buffer gas atoms can only be achieved for sufficiently low 
atom-to-ion mass ratios ^ = ma/rui. For larger mass ra¬ 
tios, the ion experiences an effective energy gain through 
elastic collisions, even if the buffer gas is at zero temper¬ 
ature, finally resulting in loss from the trap. These two 
regimes are separated by the critical mass ratio ^crit, as 
first introduced by Major and Dehmelt [M- la recent 
years, several groups have determined the critical mass 
ratio numerically [M [24] as well as analytically [26] and 
found values slightly larger than Major and Dehmelt’s 
first prediction ^crit = 1- 

In order to extend the regime of efficient cooling to 
lower Hnal temperatures and larger mass ratios, two dif¬ 


ferent approaches have been explored. Spatial confine¬ 
ment of the buffer gas, e.g. atoms stored in optical traps 
[THU], restricts collisions to the trap center where the 
micromotion is smallest, thus reducing the collision in¬ 
duced energy transfer to the macromotion El [M EH- 
Alternatively, in RF traps with higher pole orders the 
micromotion is reduced over a larger volume thus allow¬ 
ing for efficient cooling with buffer gas EH Eg. Despite 
a growing number of experiments using these approaches 
a general theoretical framework describing the influence 
of both a spatially confined buffer gas and high order RF 
traps is still lacking. 

In this paper we present a comprehensive model for 
the dynamics of a single ion interacting with a spatially 
confined buffer gas inside an RF trap of arbitrary pole 
order. The collisional kinematics can be favorably de¬ 
scribed in a reference frame assigning the micromotion 
to the buffer gas rather than the ion. Depending on the 
mass ratio the ion’s final energy is either determined by 
the buffer gas temperature (^ Ccrit) or by the effective 
energy of the ion’s micromotion (^ ^ ^crit)- For a ho¬ 
mogeneous buffer gas, the micromotion restricts cooling 
to ^ < ^crit in agreement with previous work. However, 
for a spatially confined cooling agent the emergence of an 
additional stable regime is found, thus enabling efficient 
cooling of the ion motion beyond ^crit- In this regime 
the ion’s energy distribution is determined by the energy 
of the trap’s ponderomotive potential t4ff averaged over 
the buffer gas distribution. We provide semi-analytic ex¬ 
pressions for the ion’s energy distribution for arbitrary 
mass ratios and trap multipole orders. As the averaged 
ponderomotive potential can be controlled, e.g. through 
adjustment of the trap parameters or the atom cloud size, 
the ion’s final temperature can actively be changed offer¬ 
ing perspectives for enhanced cooling {forced sympathetic 
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cooling). 

Consider an ion stored in a cylindrically symmetric RF 
ion trap of pole order n, undergoing elastic collisions with 
a spatially confined neutral buffer gas as schematically 
depicted in Fig. The ion’s motion can be separated 
into it’s micromotion and macromotion, characterized by 
the velocities Vi and Ui, respectively. The micromotion is 
an implicit function of the ion’s radial position r and the 
RF phase 


sion which is assumed to be short compared to the other 
two time scales. This holds, as long as the collision en¬ 
ergy exceeds the energy scale set by the interplay between 
micromotion and short-range atom-ion interaction, as in¬ 
troduced in [5^ , which is typically on the order of tens or 
hundreds of microkelvin. In this case, the micromotion 
remains unchanged during the collision (Fi « F), and 
the macromotion after a single elastic collision is given 
by (see Supplemental Material): 


|Fi(r)| = fc—r” ^cos$RF, (1) 

irii 

with n being the trap’s multipole order, mi denoting 
the ion’s mass and k = (7Ve“)Vb/(<j-'^o) being a con¬ 
stant depending on the system parameters: Ne~ (ion 
charge), Vg (RF-voltage), oj (RF-frequency) and Rq (in¬ 
ner trap radius). In this dynamically driven system the 
ion’s energy oscillates with $rf- Generally, this os¬ 
cillation is much faster than the time scale associated 
with the macromotion. For a Paul trap, this corre¬ 
sponds to a small g-parameter in the Matthieu equa¬ 
tions [IH]. By averaging over one RF cycle, the ion’s 
energy associated with the macromotion is obtained as 
Ei = miuf (2 -I- Ves(r) with the ponderomotive potential 
Ves = TOi(Fi(r)2)RF/2 = n2pr2"-2/4mi. 

A third time scale is defined by the duration of a colli¬ 


linear octupole trap (n=4) collision at p=0 



FIG. 1: Schematic of an ion trajectory (blue) in a linear 
octupole trap (pole order n = 4) colliding with an buffer gas 
atom (red). The velocities Vi and Ui are associated with the 
ion’s micro- and macromotion, respectively, while the atom’s 
velocity is indicated by Fa. In the lower graph, the effective 
ponderomotive potential I4fi) of the ion trap is shown to¬ 
gether with the spatial distribution of the buffer gas atoms 
(spatial extension 2aa). In the right panels, typical ion tra¬ 
jectories are shown for a collision near the the trap center, 
yielding a decrease of the ion’s energy (upper panel), and 
near the classical turning point, resulting in an increase of 
the ion’s energy (lower panel). 


A \ f- -M , Fi±i(^a-Fi) 

^ Yc) * (^a “1“ ' “ ^ j 


( 2 ) 

with Fa being the atom’s velocity and TZ{9c,(l)c) being 
the rotation matrix defined by the polar and azimuthal 
scattering angles 9c and 

Note, that Eq. ([^ is formally equivalent to an elas¬ 
tic collision in free space of an ion with momentum miUi 
colliding with an atom of momentum ma(Fa —Fi). There¬ 
fore, it is useful to choose a reference frame in which the 
micromotion is assigned to the atom instead of the ion, 
leading to an effective atom velocity Fes = Fa —Fi. By av¬ 
eraging over all collision angles and one period of the mi¬ 
cromotion the atom’s average kinetic energy is obtained 
as 


(Aa(r)) = ^ma (v^s) = ^Ves{r) + (3) 


with /cr being the Boltzmann constant and Ta the buffer 
gas temperature. 

By comparing the radial dependence of the atoms av¬ 
erage kinetic energy to the one of the ion, as given by 
^miuf{r), two distinct regions, separated by a radius Tc, 
can be identified: For r < Vc the ion’s kinetic energy ex¬ 
ceeds the average energy of the atoms resulting in a net 
energy transfer to the atoms, whereas for r > Tc the ion’s 
energy is generally increased through a collision with the 
buffer gas atom. The radius Tc is given by 


— ^max 


I 1+e J 


(4) 


with Tmax being the ion’s maximum turning point in 
the ponderomotive potential, as defined by the condition 
Veff(r’max) = E^. For E; < k^T^,, the net energy trans¬ 
fer is always positive and, thus, the radius Tc no longer 
defined. 

In order to numerically determine the ion’s equilib¬ 
rium energy distribution, the energy E^ is tracked over 
the course of many collisions. Generally, this would re¬ 
quire to solve the equations of motions for the ion’s full 
trajectory and to evaluate the scattering probability at 
every infinitesimal time step. In order to circumvent 
this computationally demanding approach, two simpli¬ 
fications can be made. Firstly, for the energy regime 
discussed here, the classical Langevin model [321131] can 
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be used, which results in a velocity-independent scatter¬ 
ing probability [24l[26]. Secondly, assuming that the time 
between consecutive collisions is long compared to the pe¬ 
riod of the macromotion, the scattering probability can 
be averaged over the macromotion resulting in an radi¬ 
ally dependent scattering probability (see Supplemental 
Material). 

Fig. i shows the resulting final distributions P{E\) 
for different mass ratios in a Paul trap. Depending on 
the mass ratio we find three different regimes char¬ 
acterized by the shape of the energy distributions. For 
^ <C 1 the ion thermalizes to the buffer gas temperature 
as the atoms’ effective energy is dominated by the atoms 
thermal energy (see Eq. @)- Therefore, the ion’s ther¬ 
mal distribution follows a Maxwell-Boltzmann distribu¬ 
tion for a gas in a harmonic radial potential (Boltzmann 
regime, see Table |l]). 

A second regime emerges for larger mass ratios ^ « 1. 
The energy distribution exhibits a power law tail E'^ to¬ 
wards higher energies, while the low energy behavior still 
follows a Boltzmann distribution {Power-Law regime, see 
Table|l| . Such power-law behavior has been identified and 
explained in previous investigations [22H2^[26] . From the 
frame transformation , resulting in Eq. (|^, the deviation 
from a Boltzmann distributions can be understood as a 
consequence of the contribution of the micromotion to 
the atoms’ effective velocity. Through collisions with the 
buffer gas, an ion can gain multiples of it’s current energy 



FIG. 2: Normalized equilibrium energy distributions for dif¬ 
ferent mass ratios ^ in a Paul trap (pole order n = 2). The 
buffer gas cloud distribution is given by a Gaussian of size 
Ca = i?o/100 and temperature of Ta = 200pK. Also shown is 
the energy distribution in the Boltzmann regime (red curve) 
and the energy distribution for ^ = 34 (purple curve), ac¬ 
cording to the expressions in Table [I] The inset compares 
the exponents k of the power-law in the energy distribution 
(as defined in Table , for different models. The condition 
K = —2 separates the regimes of stable from unstable ion mo¬ 
tion. Results found by Zipkes et al. have been corrected as 
Kcorr = K — 1 (see Supplemental Material). 


heff (t’max). 

The inset of Eig. depicts the power-law exponent 
K obtained from our simulation and compares it with 
previous models, showing excellent agreement. For n > 
—2 the ion’s mean energy diverges and the ion is no longer 
confined by the trap. This condition is commonly used 
to define the critical mass ratio ^crit [23 US] ■ From the 
simulations we obtain ^crit ~ 1.4. 

If the buffer gas is evenly distributed over the entire 
ion trap, sympathetic cooling is only feasible for ^ < ^crit. 
For an atom cloud with finite size CTa, however, we find an 
additional stable regime even for ^ ^ 1, characterized by 
an energy distribution following a power-law at lower en¬ 
ergies bound by an exponential decrease towards higher 
energies {Localization regime, see Table [I]). For increas¬ 
ing mass ratios, the exponent in the Power-Law regime 
increases, until the distribution becomes essentially flat 
{k ~ 0) for ^ « 5. For even larger mass ratios, the expo¬ 
nent increases further, finally converging towards At « 3/2 
for ^ 1. 

The exponential decrease of the energy distribution at 
the higher energies is caused by the localization of the 
buffer gas. As collisions are restricted to the volume of 
the buffer gas, the atoms’ effective energy (£'a(r)) (see 
Eq. is bound by the finite size of the cloud. The 
effective energy content Ea accessible for collisions with 
an ion is given by integrating {Eg_{r)) over the spatial 
distribution of the buffer gas. Eor a Gaussian shaped 
cloud with standard deviation Ca this yields 

^a= ^fcBTa + 2"-i(n-l)!eEeff(aa) . (5) 

In the Localization regime, the scale of the energy dis¬ 
tribution is thus given by Eg,, in contrast to k^Ta for 
the Boltzmann regime. This is illustrated in Fig. For 
^ = 0.1 the numerical data is well reproduced by a Boltz¬ 
mann distribution with energy scale k^sTa, whereas the 
distribution shown for ^ = 34 is given by a distribu¬ 
tions function of Boltzmann type, see Table [Tj {Localiza¬ 
tion regime). The mean energy is given by Eg, therefore 
the characteristic energy scale E* becomes dependent on 
the power-law exponent n. 

So far, we have discussed the case of a Paul trap. 
For higher trap orders, we find the same regimes, yet 
a critical mass ratio '^crit(?T-) depending on the multipole- 


Boltzmann regime (§ <C ^crit) P{Ei) oc E\P exp {— ) 

Power-Law regime (^ ~ ^crit) P{Ei) « | ^ 

Localization regime (^ ^ ^crit) P{Ei) oc E^ exp(— ;pr) 


TABLE I: Analytical expressions for the ion’s energy distribu¬ 
tion in the three different regimes for a Paul trap, as obtained 
from fitting the numerical results shown in Fig. 
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order n. The energy distributions are slightly modified to 
the ones described in Table |l} In the Boltzmann regime 
the distributions differ as the average energy stored in 
the ponderomotive motion varies with n. The ratio be¬ 
tween the ion’s average potential and kinetic energy can 
be expressed as (I4ff)/(^^kin) = 2/(3n — 3), as follows 
from the virial theorem. Thus, the ion has a larger 
mean energy {Ei) = (| -I- ;^)fcBTa than the buffer 
gas atoms thermal energy Therefore, the cor¬ 

responding Boltzmann distribution is given by P{Ei) oc 
^i/2-i-i/(’^-i) Qxp{—Ei/kBTa), which reproduces our nu¬ 
merical simulations well. 

With increasing mass ratio we observe the emergence 
of the Power-Law regime as in case of the Paul trap, with 
the difference that the transition now occurs at higher 
mass ratios. For a fixed mass ratio one finds that the 
exponent decreases with increasing pole order, as was 
already described in [^. The critical mass ratio, again 
defined by the exponent k = — 2, is numerically found to 
be 

^crit(»2) « 1.4(n - 1) . (6) 

The critical mass ratio can also be approximated ana¬ 
lytically from Eq. ([^ by solving {u(^ — Ui^)RF = 0 and 
again applying the virial theorem. This results in a crit¬ 
ical mass ratio of 'fcrit(fi) = 1-5 {n — 1) which reproduces 
the linear dependence on the pole order of the numerical 
findings. Consequently, Eq. suggests that sympathetic 
cooling can be efficiently applied even with a homoge¬ 
neously distributed heavy buffer gases (C ^ 1) as long as 
the pole order n is sufficiently high. 

Eor a localized buffer gas with a mass ratio exceeding 
the critical value as given in Eq. the ion’s energy dis¬ 
tribution follows the same analytic form as in a Paul trap 
(see Localization regime in Table |T]) . The corresponding 
energy scale is given by the general expression Eq. [^yield¬ 
ing P{E\) cx Ef exp(— "“-yEi/E*). The mean energy is 
again given by Ea, therefore E* depends on n as well as 
the power-law exponent k. 

The scaling of the maximum energy given by Eq. [^im¬ 
plies that the ion’s mean energy can be decreased by low¬ 
ering I4ff(o’a) {Forced Sympathetic Cooling, ESC). This 
can be achieved by either lowering the rf voltage Vq, or 
by compressing the buffer gas cloud, or both at the same 
time. Unlike forced evaporative cooling in atomic traps, 
ESC does not lead to a loss of particles [35]. Adjusting 
these parameters sufficiently slow compared to the colli- 
sional equilibration time, the ion remains trapped. We 
find in our simulations, that the ratio between the size of 
the buffer gas cloud and the volume probed by the ion re¬ 
mains constant, leaving the relative overlap and thus the 
average scattering rate unaffected. This holds, as long 
as the thermal energy in Eq. is negligible, resulting in 
the ion’s mean energy being proportional to I4ff(o’a)- In 
this case the volume probed by the ion, characterized by 


?'max, is directly proportional to tXa and does not depend 
on the trap parameters. 

Fig.i demonstrates the principle of ESC by slowly 
decreasing the size of the buffer gas. Following Eq. the 
ion’s energy distribution is shifted towards lower energies. 
In case of a Paul trap and a mass ratio of ^ = 10, this 
leads to a decrease in the ion’s temperature from 1600K 
for (Ta = 0.64i?o down to the buffer gas temperature of 
4K for tTa = 0.01i?o (see inset of Eig.[^. 

Using the favorable frame transformation, where the 
micromotion is assigned to the buffer gas atoms, we have 
found simple analytic expressions describing the energy 
distribution of the ion and provided intuitive insight into 
the dynamics of buffer-gas cooling in ion traps of arbi¬ 
trary multipole order. For a spatially confined buffer 
gas, we discovered a regime of stable ion motion even 
for large atom-to-ion mass ratios, thus extending previ¬ 
ous investigations on the dynamics of ions colliding with 
a homogeneously distributed buffer gas. These findings 
enable the design of a new class of experiments for creat¬ 
ing cold and ultracold samples of atomic and molecular 
ions. By ramping down the ion trap voltage and/or the 
spatial extension of the buffer gas, one can further cool 
the ion’s motion (Forced Sympathetic Cooling). As a re¬ 
alistic example [36] we estimate a final temperature of a 
few tens of millikelvin for OH~ ions in an octupole trap 
cooled by Rb atoms in a magneto-optical trap despite 
the mass ratio of ^ = 5. So far, we have ignored ion-ion 
interactions m as well as the axial confinement [291138] 
or the effect of stray electric fields [24]. To a first ap¬ 
proximation these additional factors might be treated as 



ion energy £, [meV] 


FIG. 3: Forced Sympathetic Cooling. Shown are four en¬ 
ergy distributions in a Paul for different buffer gas sizes. The 
buffer gas has a temperature of Ta = 4K and ^ = 10. The 
two curves correspond to the analytic expressions found in 
Table ||| for a a = 0.01 and 0.64. The inset shows the ions 
final temperature as a function of the buffer gas size. The 
analytic expression corresponds to Eq. (j^, the numeric result 
was obtained by fitting a modified Boltzmann distribution to 
the numeric data. 
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perturbations of the ponderomotive potential, resulting 
in a reduction of the critical mass ratio. 
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APPENDIX A: COLLISIONS IN AN RF TRAP - 
DERIVATION OF EQ. (2) 


Consider a single elastic collision between an atom and 
ion inside an rf trap with ma/wi and Ua/cl being their 
mass and initial velocity in the lab frame, respectively. 
All initial lab frame velocities can be expressed by the 
center-of-mass (COM) V = (miCi + m^Vg^)/(mi + ma) and 
the relative velocity vq = Ci — Va.' (ci = V + ^vo,Va = 

V + :^vo). Here /r = mima/(mi + rria) is the re¬ 
duced mass. Elastic collisions rotate the relative veloc¬ 
ity in the COM frame, without changing its’ magnitude: 
Vq = 'R{9c, 4‘c)vo, with 9q and (pQ being the polar and the 
azimuthal scattering angle, respectively. The final veloc¬ 
ity of the ion in the lab frame is obtained by adding the 
COM motion V onto the rotated relative velocity: 


rria 

m -I- rria 


TZ{9c,(l)c){ci 


Va) + 


TOjCj -I- maVa 
rrii + nia 


(SI) 


In the rf trap the ion’s motion can be decomposed into 
the macro- and the micromotion: c) = Ui+Vi. The micro¬ 
motion is a well-defined velocity only depending on the 
radial distance to the trap center r and the rf phase d)RF- 
As the duration of the collision is typically small com¬ 
pared to the period of the micromotion, one assumes that 
the micromotion remains unchanged during the course of 
the collision, resulting in c( = ui -I- Ui. The ion’s energy 
is given by Ei = ^rriiuf + VeS- The energy difference is 
given by: AE^ = — u\). The final macromotion 

velocity u{ is obtained by ci — Vi which leads to: 


-./ ^ . mUi+rUaiVa- Vi) 

Ui = -^- n[ 9 c,(pc){Ui-{Va-Vi)} + 

rrii + vria mi + ma 

jS2) 

One should note that by choosing c) = Ui and Va = 
Va — Vi, Eq. (SI) can be converted into Eq. (S2|. Conse¬ 
quently, an atom-ion collision in the driven field of an rf 
trap can be described as a collision in free space of an ion 
and atom with modified momenta miUi and ma{va — Vi). 


APPENDIX B: NUMERICAL SIMULATION OF 
THE ENERGY DISTRIBUTION 

In order to numerically calculate the ion’s energy dis¬ 
tribution, the ion is tracked over the course of typically 
10® collisions. This is commonly done by calculating the 
ion’s exact trajectory and evaluating the collision prob¬ 
ability P{t)dt for every infinitesimal time step dt. This 
approach leads to very long computation times, as the 


time between two collisions is usually much longer than 
the time scales set by micro- and macromotion. In or¬ 
der to circumvent this problem, the large difference of 
the time scales can be used to separate them. As long as 
the ion undergoes many oscillations in the ponderomotive 
potential (macromotion) between every two consecutive 
collisions, the ion’s exact trajectory can be substituted 
by a radial density distribution Pi{r). The density dis¬ 
tribution is independent of the other two cylindrical co¬ 
ordinates z and (j), as we assume perfect symmetry in 
these dimensions. In this case, the probability of an col¬ 
lision to occur in the interval dr around a radial position 
r is given by P{r) dr = Pa(v)cfv-ceiPi(r) dr, with Pa{r) be¬ 
ing the atoms’ density distribution, cr being the collision 
cross section and Urei = |wi —Wff| being the relative atom- 
ion velocity. Further assuming, that the cross section is 
given by the Langevin expression cr ex l/z^reh the relative 
velocity drops out of the scattering probability and P{r) 
is proportional to the overlap between the atoms’ density 
distribution and ion’s spatial probability distribution. 

The ion’s spatial probability distribution Pi{r) is char¬ 
acterized by the ion’s total energy and angular mo¬ 
mentum. The total energy can be separated into two 
components, the energy Elz = miuf^l2 in axial direc¬ 
tion without confining potential and the radial energy 


time [27r(o] /7,(r) 

0 100 200 300 400 1 3 5 



time [l/ioi] pi{i) 


FIG. 4: Ion’s radial probability distribution for different an¬ 
gular momenta. All angular momenta are expressed in terms 
of the the maximum angular momentum. The case L = 0 
corresponds to a trajectory through the trap center, the case 
L = 1 corresponds to a circular orbit leading to r{t) = const. 
The left graphs show the ion’s trajectory r{t), the right graphs 
show the probability distribution Pi{r) to find the ion in an 
interval dr around r. The upper graphs correspond to n=2, 
the lower graphs to n = 22. 
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Ej: = miU?j.f2 + Ves(r) consisting of macromotion and 
ponderomotive potential, with rti^r and being the ra¬ 
dial and axial components of the macromotion. Together 
with the ion’s angular momentum L = TOi|r x Ui.r|, these 
quantities are constants of the ion’s motion and uniquely 
define Pi(r). The probability distribution is proportional 
to the inverse of the ion’s radial velocity r, which can be 
expressed as 

Pi(r) oc (— {E^-Ves{r)) -. (S3) 

\mi mfr^J 

Fig .[^ shows the resulting radial probability distribution 
Pi{r) for different angular momenta. 

In the simulation, the ion’s radial and axial energy as 
well as it’s angular momentum are computed and stored 
after every collision. As initial conditions we typically 
use Er = k-QTa,, Ez = fcBTa/2 and L = 0. For every 
collision we then pick the following set of parameters: 

• RF-phase - the phase of the rf-field is ran¬ 
domly chosen between 0 and 27r. 

• Atom velocity - all three cartesian coordinates of 
the atom’s velocity vector are chosen with a normal 
distribution with standard deviation 

• Collision radius - the radial position of the col¬ 
lision Tcoii is chosen according to P(r). For large 
radial ion energies P(r) can have a long tail of near 
zero values, in which case we use an upper bound¬ 
ary of Tcoll < SCTa- 

• Scattering angles - using the Langevin model, 
the scattering angle in the center-of-mass frame is 
distributed isotropically. This is achieved by ran¬ 
domly choosing an azimuthal angle (/)c between zero 
and 27r and a polar angle dc between zero and tt tak¬ 
ing into account the Jakobian determinant, yielding 
a probability distribution p(0c) = sin 0c- 


Based on these parameters we can calculate the ions 
micro- and macromotion velocities. Plugging all of these 
parameters into Eq. (S2|, the ion’s macromotion after 
the collision is obtained which defines the ion’s radial 
and axial energy as well as it’s angular momentum after 
the collision. As a last step we compute the average time 
of free motion at this energy and angular momentum 
shell as given by the inverse of the overlap of atom and 
ion distribution t = (f Pa(r)pi(r)dr)~^. After perform¬ 
ing 10® collisions all energy values are weighted with the 
corresponding r and binned, resulting in the final energy 
distribution. 


APPENDIX C: SCALING THE POWER-LAWS 
FROM ZIPKES ET AL. [H] 


Zipkes et al. [23] determine the steady state solution 
of a single ion immersed in an ultracold atomic cloud. 
In their model the ion undergoes 10® consecutive col¬ 
lisions which are recorded and afterwards binned on a 
logarithmic scale, i.e. choosing bins with a logarithmic 
distribution of sizes. In this way, the logarithmic energy 
probability distribution dP{E\)/d\og{E{) is not properly 
normalized. The logarithmic distribution function relates 
to the normalized linear distribution dP{Ei)/dEi by 


dP{Ei)/d\og{Ei) oc El ■ dP{Ei)/dEi. (S4) 


If this operation is applied to a power-law distribution, 
the resulting power-law coefficient of the normalized dis¬ 
tribution is reduced by unity. Therefore, the coefficient 
of the power law k found in Ref. [24] has been corrected 
as Kcore = K — 1 resulting in excellent agreement with the 
other models, including ours. 





